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ABSTRACT 

In this Letter, we investigate the stabiHty of the statistical equilibrium of spherically symmet- 
ric collisionless self-gravitating systems. By calculating the second variation of the entropy, 
we find that perturbations of the relevant physical quantities should be classified as long- and 
short-range perturbations, which correspond to the long- and short-range relaxation mecha- 
nisms, respectively. We show that the statistical equilibrium states of self-gravitating systems 
are neither maximum nor minimum, but complex saddle-point entropy states, and hence differ 
greatly from the case of ideal gas. Violent relaxation should be divided into two phases. The 
first phase is the entropy-production phase, while the second phase is the entropy-decreasing 
phase. We speculate that the second-phase violent relaxation may just be the long-wave Lan- 
dau damping, which would work together with short-range relaxations to keep the system 
equilibrated around the saddle-point entropy states. 

Key words: methods: analytical - galaxies: kinematics and dynamics - cosmology: theory - 
dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

It has been realized that the conventional methods of statistical me- 
chanics of short-range interaction systems cannot be directly car- 
ried out to study long-range self-gravitating systems. Hence, it is 
necessary to return to the starting point of statistical mechanics 
and develop special techniqu es to handle the long-ran ge nature of 
IPadmanabh anl 1 9901) . In our earlier two works (IHe & Kand 
iKang & H9i201lh . we performed preliminary investigations 
of statistical mechanics of collisionless self-gravitating systems and 
found that (1) both the concept of entropy and entropy principle 
are still valid for self-gravitating systems; (2) entropy is additive 
and hence extensive; (3) the equilibrium states also consist of me- 
chanical equilibria; and (4) systems' local and global equilibria are 

different. 

Based on these findings, in iHe & Kand l201ll) . we formu- 
late a systematic theoretical framework of the statistical mechan- 
ics of spherically symmetric collisionless self-gravitating systems, 
with innovative approaches that differ significantly from the con- 
ventional statistical mechanics of short-range interaction systems. 
First, we demonstrate that the equilibrium states of self-gravitating 
systems consist of both mechanical and statistical equilibria, the 
former being characterized by a series of velocity-moment equa- 
tions and the latter by the statistical equilibrium equations, which 
should be derived from entropy principle. The velocity-moment 
equations for spherical systems of all orders are derived from the 
steady-state collisionless Boltzmann equation. Then, we point out 
that ergodicity is invalid globally for self-gravitating systems, but 
can still be reestablished locally if gravitating particles are treated 
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as indistinguishable. Based on the local ergodicity, by using Fermi- 
Dirac statistics, with the weakly-degenerate condition and the spa- 
tial independency of the microstates, we re-derive the Boltzmann- 
Gibbs entropy, which should be exactly the correct entropy form 
for collisionless self-gravitating systems. Apart from the usual con- 
straints of given mass and energy, we demonstrate that the series 
of moment equations must be included as additional constraints 
on the system's entropy functional when performing the varia- 
tio nal calculus, whic h is an extension to the original prescription 
bv lWhite & NaravanI ( fl987l) . 

In that work, we left behind a question whether the statistical 
equilibria of gravitating systems, worked out with our approach, are 
stable or not. In this Letter, we investigate the stability of the sta- 
tistical equilibria of spherically symmetric systems by calculating 
second variation of entropy. The Letter is organized as follows. In 
Section[2l we prese nt solutions with a truncated distribution func- 
tion (DF), derived in He & Kang i201 1). In Section[3] we calculate 
second variations of both ideal gas and gravitating systems to make 
a close comparison between these two systems. We present discus- 
sions and conclusions in Section|4] 



2 SOLUTIONS BASED ON THE FIRST VARIATION 



We br iefly summarize the basic equations derived in iHe & Kand 
ilOl ih . on the basis of the first variation of entropy, with the trun- 
cated DF up to the second order of velocity expansion. 
The coarse-grained DF can be formally expressed as 



F{X,V) ^ exp ( - ^ \k,m,n{x)ViV^V3), 

k,m,n 



(1) 
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where Xk^m.nix) correspond to the formal Taylor expansion coef- 
ficients of In F{x, v) with respect to v. These undetermined coeffi- 
cients should be determined by the joint mechanical and statistical 
equilibrium equations. At this moment, we do not know how to 
evaluate the entropy and derive the equations of statistical equilib- 
ria with this complete DF and hence some truncation to the DF is 
necessary. We consider that, with this truncation, although the re- 
sults are not accurate, they are still valuable and heuristic. 

For the spherically symmetric systems, the truncated DF 
F(x, v) of equation (TJ to the second order can be expressed as 



F{r,Vr,vg,v^) = 



(27r) 2 arO-ga^, 



expl 



2^ 



4> 



),(2) 



in which the density p, the velocity dispersions ct,^, and cr^ are 
all functions of r. As a result, the total (Boltzmann-Gibbs) entropy 
of the system with this truncated DF is 



5'bg = - / F In FA^cYv = 



In 



(^) 



+ -(l + ln27r) 



dr 



, 1/2 1/2 1/2 

^ 2 , iVr Pe N , 



(3) 



where in the last equality we drop the unimportant constant term, 
|(1 + ln27r), andp,, = pal, pe = pal andp^ = pcr^. 

We define the kinetic energies that are contained in the r- 
sphere of the system in three orthogonal directions as 



kr{r) 



'i'Kr'^pe{r')Ar' , 
A-Kr'^Pr{r')Ar' . 



Thus, the total kinetic energy contained in the sphere is 

kt{r) = keir) + k^r) + kr{r) 
and the mass function m{r) is 



m(r) 



47rr'^p(r')dr', 



(4) 



(5) 



(6) 



and u{r) is the potential energy of a spherical gravitating system: 
f Gpy)p(r") 3 , 3 „ 

(7) 



u{r) 



— —AttG / m{r')p{r')r'dr' , 
Jo 

where the integral is restricted to the volume of the r-sphere, Qr- 
With these variable definitions, the total mass AIt, total kinetic en- 
ergy Ek, and total potential energy By of the system are, respec- 
tively. 



Mt = m(r)\r 



Ek = fct(r)|r^oo and Ev = u{r)\r 



(8) 



From equations Q and (|6), we perform the variable transformation 
by calculating the first derivatives as 



kf, 



in which the superscript 'prime' denotes the first derivative with 
respect to r. With these ready, the total entropy of the system from 



P<l> 



Pr 



(9) 



equation ([3} is 



St 



Hi(r,m, m', kt,k[, kg, kg, K, k[.)Ar 

,,1/2,11/2,,, _ > / _ I,Ml/2 



m' In 
-I- m'ln (2 



,/5/2 



7/2 2 

7rr 



)]Ar, 



(10) 



where H\ denotes the in tegrand, the t erms e nclosed in the braces. 

As demonstrated in ' He & Kand fzOlTI) . the equilibrium states 
of self-gravitating systems consist of both mechanical and sta- 
tistical equilibria, the former being characterized by a series of 
velocity-moment equations. They are derived from the steady-state 
collisionless Boltzmann equation and should also be included as 
additional constraints on the entropy functional when performing 
the variational calculus, besides the usual constraints of mass and 
energy conservation. The second-order moment equation is just the 
familiar Jeans equation: 



dr r 



d , 28 ^ d<E> ,,,, 

dr r dr 



where vt denotes either vg or u^, and the barred quantities indi- 
cate the corresponding velocity moments. P is the usual velocity 
anisotropy parameter, defined as /3 = 1 — /vr- 

The Jeans e quation Jilt is equivalent to the fol lowing virial- 
ization form (e.g. lMo. van den Bosch & Whitell2010l . p.235): 

2fct(r) + u(r-) - 2rfe;(r) = 0. (12) 

We will use this relation, instead of the original Jeans equation Jilt , 
as the constraint of the mechanical equilibrium. Additionally, from 
equation Q, we have 

Gmm' 



-AnGmpr - 



(13) 



where the superscript 'prime' again indicates the first derivative 
with respect to the radius r. The mathematical relation between 
u and m of equation l ll3t provides another additional constraint. 

With all these constraints considered, the total constrained en- 
tropy of the system is 



5*1, c — 



H2{r, m, m', kt, k[, kg, kg, kr, fcj., u, u')Ar 

kg kr (^k^ kg /Cy. ) ^ 







LI 


m' In 







1/5/2 



+ m' ln(2^/^7rr^) + /i (r)(2fet +u- 2rk'r) 
+ /2(r)(u' + — — ) Ur, 



(14) 



where H2 denotes the integrand enclosed in the braces, and /i and 
/2 are two Lagrangian multipliers. By using the variable transfor- 
mations of equations l|5j and l|6), the mass and energy conservation 
constraints for the variational calculus are converted to satisfying 
the fixed end-point conditions of equation For spherical sys- 
tems, pg{r) — P4,{r), so kg{r) = k^{r). Then perform the stan- 
dard variational calculus, 5St,c = 0, with respect to the variable 
pairs (m, m'), {kt, fcQ, {kr, k'r), and {u, it'), we obtain the follow- 
ing two equations: 

dlnp _ dlnpt _ p^ _|_ 1 _ g 
dr dr rpr r ' 



3 d In p 1 Apr 2p 



2 dr 



Pt dr rpt ^ r 



Gm 



(15) 
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where pt denotes either pe or ; 
phers as 



/i 



2rpr 



P 

2rpt ' 



/2 = 



, with the two Lagrangian multi- 
P 



2pt 



A, 



(16) 



whe re A is an integration constant and s hould be always negative 
(see iHe & Kandl201Cl : [Kang & Hdl20lt . 



3 STABILITY BASED ON THE SECOND VARIATION 

The above equations are just the lowest-order approximation, de- 
rived from the truncated DF of equation (|2), but, as mentioned pre- 
viously, we may still acquire useful knowledge with this approxi- 
mated treatment. We investigate the stability of statistical equilibria 
of the self-gravitating systems by analysing the second variations of 
entropy. 



3.1 Ideal gas 

We put an ideal gas in a spherical container and perform the varia- 
tional calculus in a similar way to that for spherical self-gravitating 
systems to make a close comparison between the two systems. 

Since the Boltzmann-Gibbs entropy is suitable for short-range 
interaction systems, so the entropy of equation dlO) can also be ap- 
plied to ideal gas. We perform the first variational calculus with 
respect to all the variable pairs (m,m'), {kt,k[), {kg, kg) and 
{kr, k'j.), and obtain p — po, and Pe ~ P4. = Pr = Po, in which 
po and Po are two constants. We then calculate the second variation 
and obtain 



1 / 

^''S't = - / [Ai(Smf + Bi{Sm'f + A2{5ktf + B2{5ki) 
^ Jo 



/^2 



+ A3{SkgY + BsiSkeY + Ai{Skrf + B4iSk'.Y]dr, (17) 



where 



A 



(18) 



dy^ dr^dyidy'/' dy'^ 

in which (yi,y'i) indicates the variable pairs (m,m'), {ki,k[), 
{kg, kg) or {kr,k'^). For ideal gas, all the yl-terms are vanishing 
and the _B-terms are 

B3 _ B4 _ po 



Bi 



<0, B2 



< 0. (19) 



SvrrVo " " 2 2 2-Kr'^pl 

These results are listed in Table [T] 

We point out that, as shown in fig. 4 o f lHe&Kan8ll2010l) . 
and Sm' (or 5p) are two different modes of perturbations. The for- 
mer is the long-range (or large-scale) perturbation, whereas the lat- 
ter is the short-range (or small-scale) perturbation. Similarly, 5kg, 
5k^ and 5kr are all long-range perturbations of the relevant quan- 
tities, and 5k'g, 5k'^ and 5k'r (or 5pg, Sp^ and 5pr) are the corre- 
sponding short-range perturbations. The long- and short-range re- 
laxation mechanisms correspond to the long- and short-range per- 
turbations, respectively. We know that the relaxation processes for 
ideal gas are collisions between the particles and between the par- 
ticles and the walls of the container, which are both short-range 
relaxations. That is the reason why all the A-terms, representing 
the long-range relaxations, vanish. 

Since the _B-terms represent short-range relaxations, and they 
are all negative, so the system at every local volume element takes 
the local maximum entropy and the system acquires its cumulative 
maximum entropy from all these local maxima. 



Table 1. A- and _B-terms of the second variation of entropy for both ideal 
gas and self-gravitating systems. 



Variables 




Ideal gas 


Self-gravitating system 


(m, m') 


Ai 


= 0, Si < 


Ai > 0, Bi < 


(kt,K ) 


A2 


= 0, B2 <0 


A2 = 0, B2 <0 


(k0,k'g) 


A3 


= 0, B3 < 


A3 = 0, B3 <0 




A4 


= 0, B4 < 


A4 = 0, B4 < 



3.2 Equilibrated self-gravitating systems 

The constrained entropy of equation ( I14t is incomplete, since we 
just used the truncated DF of equation (|2}, and only with the 
second-order Jeans equation (in its virialization form) included as 
the additional constraint. However, we believe that the following 
results based on such a truncated DF are still heuristic. The second 
variation of the entropy of spherical self-gravitating systems is 



5^ St, 



[Ai{5m)^ + Bi{5m')^ + A2{5kt)^ + B2{5k[)^ 

+ A3{5kg)^ + B3{5k'gf + Ai{5krf + Bi{5k',f 

+ Ai{5u)'' + B,,{5u'f]dr, (20) 

where Ai and Bi for (m, m'), (fct, k[), {kg, kg) and {kr, k'^) are 
calculated in the same way as those for ideal gas in equation dlSt . 
except that Hi is replaced by H2 here. The non-vanishing terms 
are 

5 

^11 = 

„2 



Ai 

B2 = 
and 

As = 



Pt 
B3 
2 



P 

2pr 



Bi 



Snr'^p 



<0, 



2-!Yr^p^ 



<0, 



Ba = B2 — 



2iTr^Pr 



<0 (21) 



d'H2 



d , d'^H2 



= 0, 



d^H2 



= 0. 



(22) 



2 dr ^ dudu' 
From the s i mulat ion results of iNavarro et al.l fcoiol) and 
iLudlow et al we know that pt 55 Pr, and with A < 0, 

we have Ai > 0. These results are also listed in Table[T] 

We see that A\ > 0, but all the _B-terms are negative. Hence 
the statistical equilibrium states of self-gravitating system are nei- 
ther maximum nor minimum, but saddle-point entropy states. As 
analysed previously, 5m is the long-range perturbation, and A- 
terms represent long-range relaxations, whereas 5p, 5pr and Spt 
are short-range perturbations, and all the B-terms represent short- 
range relaxations. Thus, the saddle-point solution suggests that the 
equilibrium entropy is the cumulative maximum entropy from ev- 
ery local maximum under short-range relaxations, while simulta- 
neously being also the global minimum entropy under long-range 
relaxations. Such a saddle-point solution is completely different 
from the c ase of ideal gas. T his result also confirms our earlier 
finding of He & Kand <201Cl) . in which we performed a prelimi- 
nary study , by just u s ing a phenomenological entropy of ideal gas 
of I White & NaravanI 119871) . 



3.3 Free collapsing self-gravitating systems 

The above discussions based on the second variations are under the 
presumption that the system has already settled into equilibrium 
states, that is, 5St,c = 0. However, these discussions are invalid if 
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Equilibrium Entropy 



Violent Relaxation 
(second phase) 




Short-Range Relaxation 



Violent Relaxation 
(first phase) 



Figure 1. Illustration of the saddle-point entropy states of equilibrated self- 
gravitating systems and the two-phased nature of violent relaxati ons. T he 
saddle-point entropy solutions are a lso consistent with Antono vl h962h 's 
proof and Binney's argument jBinnev & TremaindbOOSl) that the equilib- 
rium states of self-gravitating systems are not maximum entropy states. 



the system has not arrived at its equilibrium states, say, is just un- 
dergoing free collapsing. In this case, the entropy variation should 
be analysed with its first variation. 

We assume that a free-falling isolated spherical system has a 
given mass Mt and energy Et- In this case, the virilization re- 
lation, equation l ll2t . should not be applied here. Thus, the total 
conserved energy is the only constraint, given by 



Et ~ Ek + Ev 



' dr[{k'o + k'^ + K)-^^], (23) 



and the total constrained entropy from equations (|3) or (1101 for the 
free collapsing system (FCS) is, 

/•oo 

Jo 



m' In 



+ \{kQ + k^ + kj. 



i.'l/2i.'l/2j,/l/2 

(p 

m'5/2 
Gmm' , 



+ m' ln(2^/^7rr^) 



)Ur, 



(24) 



where denotes all the terms enclosed in the braces, with A be- 
ing the Lagrangian multiplier, and should always be negativfl We 
calculate the first variation of the entropy as 



/•oo 

SSfcs = / dr[CiSm + DiSm' + C25ke + D2Sk'o 
Jo 

+(73<5fc<^ + DiSk'^ + CaSK + Di5k'^] , 



(25) 



where 
Ci = 

C2 = 

and 



am r 

akg ak^ 



dkr 



(26) 



^ Generally, the DF derived from statistical mechanics is always a function 
of the energy level e, as /(e) ~ e^^, in which we see that /(e) — > oo as 
e — >■ oo if A is positive. Hence, to protect the DF from 'ultraviolet catastro- 
phe', A should always be negative. 



dm' r ^ p'^'-' ' 



2' 



D2 = — = — + A Da = — ^—^\ 
dk'g 2k'g dk'^ 2k'^ 

= ^ = II}L + x (27) 

Note that in the free-collapsing process, especially in the starting 
period, the mass contained in the r-sphere increases dramatically 
and hence for any r, the long-range perturbation 5m ^ 0, but 
all the short-range perturbations, 5p, 5pg, Sp^ and Sp,., compared 
with Sm, are small and negligible. Thus, according to the above 
analysis, we have SSfcs ~ J drCiSm > and we arrive at the 
conclusion that the free collapsing of self-gravitating systems is an 
entropy-production process. The entropy production will decrease 
or may stop when Sm ~ 0. 

We put together the results of both sections 13.21 and 13.31 in 
Fi g. [T] The nom enclature of two-phased violent relaxations is af- 
ter ^oken d 19961) who had already realized that the entropy of self- 
gravitating systems ma y not necessarily increase during the vio lent 
relaxation proces^ (cf.'Tremain e. Henon & Lvnden-Belll 19861 : see 
also ,Sridhar..l987. : ,Kandrua.l98C 



4 DISCUSSIONS AND CONCLUSIONS 

In this Letter, we investigate the stability of the statistical equi- 
librium of spherically symmetric collisionless self-gravitating sys- 
tems by calculating the second variation of entropy. The entropy, in 
the Boltzmann-Gibbs form, is approximated by truncating the DF 
to the second-order expansion of velocity, with the constraints of 
the given mass, energy and the virialization relation, equation il2l . 
The main results of this Letter are as follows: 

(i) Perturbations of the relevant quantities should be divided into 
two categories: long-range perturbations, such as Sm, and short- 
range perturbations, such as Sp, Spr and Spt. 

(ii) Relaxation mechanisms of different scales should be rele- 
vant to the perturbations of the same scales. The long-range relax- 
ations, such as violent relaxation and the long-wave Landau damp- 
ing, are responsible for suppressing the long-range perturbations, 
whereas the short-range relaxations, such as phas e-mixing, chaotic 
mixin g or short-wav e Landau damping (see Binnev & T remairi3 
i2008i : lMo et al.l2010l) . may be responsible for wiping out the short- 
range perturbations. 

(iii) The statistical equilibrium states of self-gravitating sys- 
tems are neither maximum nor minimum, but saddle-p oint entropy 
states. This finding conforms to ou r earlier r e sult o f iHe & Kand 
( |2010[) and is also consistent with lAntoiiov ( Il962h 's proof and 
Binney's argument jBinnev & Tremairi3l2008h that the equilibrium 
states of self-gravitating systems are not maximum entropy states. 

(iv) Violent relaxation can be classified into tw o phases. The 
two-phased violent relaxation was first proposed by ISokeij d 19961) 
who had already realized that the entropy of self-gravitating sys- 
tems may not necessarily increase during the violent relaxation pro- 
cess. Our saddle-point solutions support this classification. 



The free collapsing of gravitating systems, driven by (the first- 
phase) violent relaxation, is a rapid entropy-production process. 



The violent relaxation we referred to in fig. 5 of 'He & Kan j i20\(jt) 
should actually be the second-phase violent relaxation of this Letter 
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When the long-range perturbation Sm vanishes, the violent relax- 
ation will cease and the rapid entropy production may also stop. 
However, the entropy might still increase slowly, driven by the 
short-range relaxations to approach the final equilibrium states, 
while if relaxations were so efficient that the produced entropy 
would exceed the equilibrium entropy, then the second-phase vi- 
olent relaxation might be induced to draw the system back to its 
equilibrium state by decreasing its entropy. Such a process, that 
is, oscillating around the equilibrium entropy, would last for some 
time with a gradually damping amplitude and finally the system 
would arrive at its equilibrium state. 

The above description, suggested by the saddle-point solu- 
tions, may be a possible picture of the stability of the statistical 
equilibrium for self-gravitating systems and we further speculate 
that the second-phase violent relaxation may just be the long-wave 
Landau damping. The justification of the picture as well as the 
speculation need further investigations. 
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